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The theoretical aspects of an orbit determination Biter that incorporates ground- 
system error sources as model parameters for use in interplanetary navigation are 
presented in this article . This filter , which is derived from sequential filtering theory, 
allows a systematic treatment of errors in calibrations of transmission media, station 
locations, and Earth orientation models associated with ground-based radio metric 
data, in addition to the modelling of the spacecraft dynamics. The discussion 
includes a mathematical description of the filter and an analytical comparison of its 
characteristics with more traditional filtering techniques used in this application. 
The analysis in this article shows that this filter has the potential to generate 
navigation products of substantially greater accuracy than more traditional filtering 
procedures. 


I. Introduction 

In JPL’s interplanetary orbit determination process, 
ground error sources associated with radio metric data, 
such as station location, Earth orientation, and transmis- 
sion media calibration errors, are usually treated as un- 
modelled “consider” parameters in a scheme known as the 
consider option. This method does not utilize any infor- 
mation pertinent to the consider parameters in computing 
estimates of the trajectory parameters. Rather, the effects 
of the consider parameters are accounted for by modify- 
ing the computed estimation error covariance with pre- 
assigned uncertainties of these parameters. In many cases, 
this method gives satisfactory results and allows reason- 
able navigational accuracy. However, the consider option 
sometimes yields unstable and unpredictable results when 
used for interplanetary navigation. 


Recently, an enhanced orbit determination scheme 
(which will be referred to throughout this article as the 
“enhanced” filter) has been developed. The enhanced fil- 
ter explicitly models ground error sources as random pro- 
cesses simultaneously with the trajectory-related parame- 
ters. This method exploits the full information content of 
the data pertaining to these ground errors in the filtering 
process and may thereby improve the knowledge of these 
parameters and the overall accuracy of the estimated flight 
path. 

The enhanced filter has been successfully applied, to se- 
lected navigation problems in some interplanetary orbit- 
determination case studies at JPL [1,2]. These studies 
have shown an increase in orbit-determination accuracy 
using the enhanced filter method of factors of two to four 
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over the use of more conventional techniques. Interest- 
ingly, in the field of Earth satellite orbit determination, 
similar filtering techniques have been in use for some time. 
Numerous articles exist in the literature. For example, in 
the article by Lichten and Border [3], stochastic models 
for the tropospheric delays are used for problems in Global 
Positioning System orbit determination. 

This article describes the framework of the enhanced 
filter model and gives an analytical comparison of the en- 
hanced filter with the traditional consider option. Since 
there are many concepts involved, a quick review of the 
basic filter model from which the enhanced filter evolved 
will be given first. A detailed discussion of the consider 
option will also be presented, so that its characteristics can 
be compared with those of the enhanced filter. 

II. Basic Filter Model 

All of the filters discussed herein employ a linear repre- 
sentation of the process dynamics and the measurements. 
In this formulation, the state space and measurements are 
described as follows: 


Xj + l = $(j + l»i)X; + Wj 

(1) 

z j — Hj Xj + Vj 

(2) 


where x denotes the extended state, which includes the 
spacecraft state and the dynamic process perturbation pa- 
rameters. The dynamic process perturbations may be 
treated as random or deterministic parameters. The tran- 
sition matrix from time tj to tj+ \ is $(j -f 1 , j), and w 
represents the state-space modelling errors. The j th data 
point is z j y and Vj is the corresponding data noise. The 
term H j is a matrix that contains the partial derivatives 
of the jth data point with respect to the extended state. 

The computations for the estimate and c ova riance ma- 
trix follow the basic square-root information filtering pro- 
cedure discussed by Bierman [4]. This procedure can be 
described as follows: The a priori values and covariance 
matrix for the estimated parameters are first transformed 
into the a priori residual and the a priori information ma- 
trix. To ensure numerical stability, the Cholesky decom- 
position is applied to the a priori information matrix to 
obtain an upper triangular factorization, hence the term 
“square root” of the matrix. This a priori square-root in- 
formation matrix is augmented by the a priori residual. 
Measurements are then included into this augmented ma- 
trix using the Householder orthogonal transformations, re- 


sulting in the following a posteriori square-root informa- 
tion and residual matrix: 



where R denotes the information matrix, z denotes the 
residual, and e denotes the sum of squares of the residual 
errors as defined in the classical least-squares problem. 

The extended state estimate x and the error covariance 
matrix for the extended state are computed using the 
values in the above matrix. 


R -1 z 

( 4 ) 

R -1 R -t 

( 5 ) 


III. Consider Option 

The consider option is based on the assumption that 
there is a set of parameters that affects the performance 
of the filter and that it is unnecessary or impractical to 
model these parameters accurately. These parameters are 
referred to as consider parameters. Since consider parame- 
ters are not included in the model, they are not estimated. 
Moreover, the estimation and covariance computation pro- 
cess for the extended state is not aware of the presence of 
these parameters and any errors in their values. Instead, 
once the state estimation is performed and the error covari- 
ance computed, the consider filter modifies the computed 
error covariance to account for a constant uncertainty in 
the consider parameters. This new covariance is called the 
“consider covariance.” If P x denotes the computed error 
covariance matrix and P yo denotes the a priori covariance 
of the consider parameters, then the consider covariance, 
P c , is computed as follows: 

P c = Px + SP !/0 S t (6) 

The matrix S above is commonly called the sensitivity ma- 
trix and contains the partial derivatives of the estimated 
state with respect to the consider parameters. In the con- 
sider filter, it is a function of the measurements only. More 
details of the sensitivity matrix can be found in [4]. 

The consider option is often used to treat ground error 
sources that affect the measurements. The rationalization 
is that since these error sources do not affect the spacecraft 
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(7) 


dynamics, estimating them will not improve the knowledge 
of the spacecraft state, and therefore it is adequate to just 
characterize them by their uncertainties. In this approach, 
the presence of the ground errors is acknowledged in a con- 
servative fashion, which overlooks the important fact that 
these error sources, since they affect the measurements, do 
affect the spacecraft state estimates. Note also that this 
method assumes that the ground error parameters are con- 
stants with no dynamics. This simplified modelling does 
not allow an improvement in knowledge of these parame- 
ters to be obtained. 


yj+i = *0' + i,i)y; + u i 

where y denotes the vector consisting of the ground error 
parameters, \P(j + 1 ,j) accounts for the time dependency 
of y from time tj to /j + i, and u represents the random 
driving term. 

The measurement model given by Eq. (2) is modified 
to include information pertaining to the ground error pa- 
rameters. 


The most obvious disadvantage of the consider option 
is that when the sensitivity of the state with respect to the 
ground error parameters increases, the consider covariance 
increases, as can be seen clearly from Eq. (6). This im- 
plies a greater uncertainty in the spacecraft state. When 
this happens, the only remedy available within the con- 
text of the consider option is to decrease the “weight” of 
the data, i.e., assume that each measurement contains less 
information than it does in actuality. This and other char- 
acteristics of the consider option are discussed by Scheeres 

[5]- 

There is, however, a deeper significance here. If the 
sensitivity of the state with respect to the ground error 
parameters grows large, then this implies that there is sig- 
nificant information contained in the measurements con- 
cerning these parameters. The logical recourse would be to 
exploit this information in order to learn something about 
the characteristics of the ground error sources. This im- 
plies that the ground error parameters should be incorpo- 
rated into the filter model. 


Zj = Hj Xj + Gj y j + Vj (8) 

where the meanings of H ; and Vj are the same as before 
and Gj contains the partial derivatives of the jth data 
point with respect to the ground error parameters. 

The enhanced filter model thus consists of Eqs. (1), (7), 
and (8). The characteristics of this model can be described 
as follows. Ground error sources are treated as system pa- 
rameters that can be estimated. The extended spacecraft 
state space evolves independently from the ground system 
parameters. The measurements have explicit dependency 
on the ground system parameters. 

B. Estimate and Covariance 

In the enhanced filter, estimated parameters include 
both the extended spacecraft state and the ground system 
parameters. The computations for the estimate and co- 
variance matrix follow the basic square-root information 
filtering procedure described in Section II. After process- 
ing measurements, the a posteriori information and the 
residual matrix have the following form: 


IV. Enhanced Filter 

The enhanced filter provides for the inclusion of the 
ground error parameters in the filter model so that the 
estimation process incorporates the information pertinent 
to the behavior of these parameters. In this treatment, the 
ground error sources are modelled as dynamic entities. By 
doing this, an automatic feedback mechanism is created 
so that no artificial data weighting is needed to keep the 
covariance of the state estimate from diverging. 

A. Enhanced Filter Model 

The enhanced filter uses the extended state model for 
the basic batch sequential model given by Eq. (1). In ad- 
dition, there are models for ground error parameters and 
a modified measurement model. The ground error param- 
eters are modelled as a discrete first-order Markov process 
in a vector form given by 


f Ibjp Rj;y 7t X \ 

H y Z y 

\ e / 


(9) 


In Eq. (9), subscripts are given to show the relationship of 
the quantities in the matrix with respect to the parame- 
ters. Here, R x represents the information with respect to 
the extended state and is the same as R in Eq. (3); R^ rep- 
resents the information with respect to the ground system 
parameters; R ry denotes the information concerning the 
extended state affected by the ground system parameters; 
and z x and z y are residual components corresponding to 
the extended state and ground system parameters, respec- 
tively. From this matrix, the extended state and ground 
system parameter estimates are computed as 
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x= R r 1 [z x - R rj , y] (10) 

y = R-y 1 Xy (11) 

Equations (10) and (11) show that the information con- 
tent pertinent to the ground system parameters, R y and 
z y , is used to estimate these parameters, yielding a new 
estimate, y. This estimate is then used to obtain a new 
state estimate, x, together with the information pertinent 
to the state, R r , the state with respect to the ground 
system parameters, R and the residual, z x . In the con- 
sider option, the state estimate would have used only R x 
and z r , and no estimate for the ground system parameters 
would have been computed. 

Since all parameters are estimated, the covariance ma- 
trix for the extended state is then 

P = RJ 1 RJ T + [~Rr *R*yRy '] [-R;‘R*yV] T (12) 

Recognizing that -K~ l H xy is the definition of the sensi- 
tivity matrix, Eq. (12) may be represented as 

T = P x + SP y S T (13) 

In this equation, P r = R-'Rj 7 and Py = R~ 'Ry 7 ”, 
where P x and P y represent the blocks of the computed 
covariance matrix corresponding to the extended state and 
the ground system parameters, respectively. 

Using the same notations, the full covariance matrix 
corresponding to the estimate for both the extended state 
and the ground system parameters is 

P^ + SPyS 7, SPy\ 

(14) 

PyS T Py / 

Note that the P x in both Eqs. (13) and (14) is ex- 
actly the computed covariance in the consider option. 
Equation (13) has the same form as Eq. (6), which gives 
the computation for the consider covariance. But the 
important difference here is that the covariance for the 
ground system parameters, P y in Eqs. (13) and (14), is 
updated using the measurements and their modelled be- 
havior, while in the consider option only the constant a 
priori covariance is used. 


Further, in the enhanced filter, the full covariance ma- 
trix is used to automatically weight the measurement data, 
thus using the modelled behavior of ground system param- 
eters to control the relative weight given to the measure- 
ments. In the consider option, only the computed covari- 
ance P x is used to weight the measurements; therefore, 
the ground system parameters have no effects on the data 
weights. In particular, using the enhanced filter formula- 
tion, the computed covariance stabilizes as the informa- 
tion contained in Ry increases. This is transparent since 
P y = R“ i R~ T . Thus, there is no need to artificially 
deweight the data when the information content of the 
ground system parameters is large. On the other hand, 
when there is little or no information pertinent to the 
ground system parameters, the enhanced filter then be- 
haves as in the consider option. In this case, the covariance 
matrix with respect to the ground system parameters, P y , 
will consist primarily of the a priori constant covariance 
for these parameters. 

C. Features of the Enhanced Filter 

The advantage of modelling measurement error sources 
as system parameters is that this makes it possible to esti- 
mate these parameters. Thus, the knowledge of them may 
be improved. This approach distinguishes the enhanced fil- 
ter strategy from other commonly used schemes that treat 
ground system errors as measurement noise. Many prob- 
lems arise with such schemes. For example, because of the 
very assumption that the ground system errors are noise- 
like, no information regarding these parameters can be ex- 
tracted from the measurements. Very often, this leads to 
degenerate covariance matrices; see [6,7] for more discus- 
sions on such treatments. 

The advantage of modelling the evolution of the ground 
system parameters independently from the state dynamics 
is that, while not affecting the state evolution, the ground 
errors are allowed to evolve according to their own dy- 
namics. This is necessary since some of the error sources, 
such as Earth orientation and transmission media, are dy- 
namic. This model demonstrates that system parameters 
do not have to be included in the state dynamics in order 
to improve the state estimation. Finally, the advantage 
of modelling the dependency of the measurements on the 
ground system parameters explicitly is that the best pos- 
sible weighting of the data can be utilized in the filtering 
process to generate estimates of the spacecraft trajectory. 

The net combination of the above features provides a 
procedure that fully exploits the information content of the 
measurements pertaining to the ground system parame- 
ters and systematically assigns the proper weight for each 
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measurement according to the modelled behavior. This 
is achieved in the enhanced filter by using user-input pa- 
rameters that describe the stochastic nature of the ground 
system parameters. By choosing the statistics of these pa- 
rameters properly, accurate estimates of both the ground 
system parameters and the spacecraft state may be ob- 
tained. In addition, the covariance of the ground system 
parameters may be effectively controlled, ensuring that the 
estimates for these parameters will be within the accuracy 
with which they can be modelled. 


V. Conclusions 

To summarize, the enhanced filter has been shown to 
offer some advantages over the traditional consider option 
in the following ways: First, the enhanced filter allows the 
behavior of the ground system parameters to affect the 
spacecraft state estimates. This helps to ensure that the 
spacecraft state is being estimated with a more complete 
representation of the physical world. Second, the enhanced 


filter can exploit the information contained in the data per- 
taining to the ground system parameters, possibly helping 
to improve the knowledge of these parameters. Moreover, 
the improved knowledge is fed back into the filtering pro- 
cess automatically, which effectively adjusts the weighting 
of the data systematically, helping to stabilize the state 
covariance. Third, the performance of the enhanced filter 
is no worse than that of the consider option in the case 
when the information content of the data with respect to 
the ground system parameters is small. 

In reviewing the history and evolution of the sequen- 
tial filtering techniques, a major motivation for using the 
consider option in the past was that it was cost effective, 
in terms of computation time, to deal with a lower dimen- 
sional model, even if the accuracy of the filtering product 
was compromised. With modern computers, this moti- 
vation is no longer a valid concern. The enhanced filter 
model can use this computational power to achieve an un- 
precedented degree of accuracy and robustness in many 
orbit determination problems. 
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